Thermodynamics of Local Causal Horizons 
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We propose an expression for the entropy density associated with the Local Causal Horizons in 
any difFeomorphism invariant theory of gravity. If the black-hole entropy of the theory satisfies the 
physical process version of the first law of thermodynamics then our proposed entropy satisfies the 
Clausius relation. Thus, our study shows that the thermodynamic nature of the spacetime horizons 
is not restricted to the black holes but it also applies to the local causal horizons in the neighborhood 
of any point in the spacetime. 



Black-hole physics has provided strong hints of a deep 
and fundamental relationship between gravitation, ther- 
modynamics and quantum theory. At the heart of this re- 
lationship is black-hole thermodynamics, which says that 
certain laws of black hole mechanics are, in fact, simply 
the ordinary laws of thermodynamics applied to a system 
containing a black hole \l,- Classical and semi-classical 
analyses of the thermodynamic behavior of black holes 
has given rise to most of our present physical insights into 
the nature of quantum phenomena occurring in strong 
gravitational fields. For example, the holographic princi- 
ple emerged through various thought experiments related 
to the sub-extensive scaling of Bekenstein-Hawking en- 
tropy - and applied it to arbitrary gravitational systems. 

The equilibrium state version of the first law of black- 
hole thermodynamics for arbitrary diffeomorphism- 
invariant theory of gravity was established by Wald and 
collaborators [21 [3] . The entropy of the black hole can be 
expressed as a local geometric quantity integrated over a 
spacelike cross section of the horizon and is identified 
with the Noether charge corresponding to the Killing 
isometry that generates the horizon. It is also possible 
to write down a quasi-stationary version of the second 
law for Lanczos-Lovelock gravity [U |S] for physical pro- 
cesses in which the horizon is perturbed by the accretion 
of positive energy matter and the black hole ultimately 
settles down to a stationary state. 

The discovery of the black-hole thermodynamics was 
also a precursor to the idea of emergent gravity, which 
asserts that spacetime is an emergent notion, a macro- 
scopic and coarse-grained description of the underlying 
degrees of freedom which are fundamentally described 
by a quantum theory of gravity. But how could gravita- 
tion be an emergent phenomenon if the thermodynamics 
is tied to special spacetimes - those containing a black 
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hole? In this direction a profound insight was provided by 
Jacobson who managed to derive the Einstein equation 
as a consequence of the thermodynamics associated to 
the Local Causal Horizons (LCHs) [6]. Jacobson showed 
that if one identifies a multiple of area with the entropy 
[S (X A) associated with a LCH at boost temperature 
T, and if the change in the entropy 6S is related to the 
flow of the boost energy 6Q across the local horizon in 
such a way that the Clausius relation dS — SQ/T is sat- 
isfied, then it is possible to derive the Einstein equation. 
The Einstein equation can therefore be interpreted as a 
thermodynamic equation of state. 

This simple and elegant result has far reaching conse- 
quences. It supports the idea that the dynamics of space 
time can be interpreted as a thermodynamic description 
of some underlying degrees of freedom. Also, as pointed 
out by Padmanabhan [7], one can deduce the presence 
of these underlying degrees of freedom by the very ob- 
servation that one can heat the system. The existence 
of physical quantities which provide a thermodynamical 
interpretation of the geometry of spacetime and the ex- 
istence of a thermodynamical law (the Clausius relation) 
that gives rise to the dynamics of the geometry (the Ein- 
stein equation) is highly suggestive of the idea that the 
dynamical geometry itself could be an emergent notion 
borne out of the dynamics of some underlying degrees of 
freedom. 

On very general grounds, however, we know that Ein- 
stein equation is not the end of the story. The Einstein- 
Hilbert action has just the first of the terms in a deriva- 
tive expansion and one expects d la Wilson that higher- 
order curvature terms respecting the symmetries should 
also appear in the effective action. It is natural to ask if 
the thermodynamical nature of LCHs is unique to general 
relativity or is it a robust feature of any diffeomorphism 
invariant theory of gravity. Although there have been 
many attempts [5l413j to answer this question and extend 
the result of Ref. [6| to other theories of gravity, none of 
these are completely satisfactory. Either they work only 
for a very specific class of theories or they require extra 



assumptions. A detailed critique and the limitations of 
these attempts is discussed in Ref. [H]. 

In this paper we take the point of view that if the 
thermodynamics of LCHs is a fundamental property of a 
theory of gravity and there are internal degrees of free- 
dom whose coarse-grained description is provided by the 
low energy gravitational effective action then we must 
be able to show the existence of thermodynamical laws 
associated to the LCHs with appropriately defined quan- 
tities as the thermodynamic variables. Hence, similar to 
black hole thermodynamics, the LCH must obey vari- 
ous thermodynamical laws. Thus, our goal is opposite to 
the one that Jacobson had in Ref. [5] . Instead of deriving 
the field equation of gravity from the assumed thermody- 
namics associated to the LCHs, we use the field equations 
to derive the Clausius relation for the LCHs. More pre- 
cisely, for any given diffeomorphism invariant theory of 
gravity in which black holes satisfy the physical process 
version of the first law of thermodynamics, we propose 
an expression of the entropy associated with the LCHs 
that satisfies the Clausius relation. In other words, we 
generalize the thermodynamics of LCH from general rel- 
ativity to any arbitrary diffeomorphism invariant theory 
of gravity. 

Geometric setup. — We start with the construction of 
the LCH based at any point p in spacetime. We first 
choose a {D — 2)-dim spacelike surface E containing p 
and we choose one side of the boundary of the past of E. 
In a small neighborhood of p this boundary is generated 
by a congruence of null geodesies fc" normal to E. Let 
A be the affine parameter along the integral curves of /c° 
such that X — Xf a.t p. Our choice of the patch E is such 
that the expansion 9 and the shear aab of fc" vanish at 
p. Our chosen side of the null boundary of the past of 
such a E is called the local causal horizon (LCH)^ at p, 
see Fig. [T] 

In order to define the heat flux and the temperature 
we need an approximate boost Killing vector field. Con- 
sider a point pq in the past of p lying at the value of the 
affine parameter A = on the central generator (i.e., the 
generator that passes through p) . Next we construct an 
approximate boost Killing vector field ^° such that the 
point Pa serves as a bifurcation point, i.e., ^° vanishes 
at Pq. On the central generator the approximate boost 
Killing vector is related to the null-geodesic generator as 
£,"' = Afc". Furthermore, ^° satisfies the Killing equation 
V(Q^b) = exactly at p, and to 0{x) near p. We suggest 
the interested reader to consult the references [131 [H] for 
further details pertinent to this construction. The appro- 
priate temperature is T = l/27r which is the boost tem- 
perature associated with the approximate Killing vector 



^ Our construction of the LCH is different from the one orginally 
given by Jacobson [6]. There the equihbrium was at the bifurca- 
tion point (where the Kilhng vector vanishes), while for us the 
equihbrium point (p) lies in the future of the bifurcation point 
(po)- Both of these constructions are discussed in Ref. I13| . 





Figure 1: Local Causal Horizon (LCH) associated with the 
approximate boost Killing vector field ^". po is the bifurca- 
tion point where ^° vanishes, fc" is the affinely parametrized 
tangent to the null generators of the LCH with the affine 
parameter A. p is the final equihbrium point where the ex- 
pansion and shear a of fc" vanishes. The affine parameter 
vanishes at po and takes a value A/ at p. On the generator 
connecting po to p we have ^^ = Afc". 



f . The heat flux is given in terms of the matter stress 
tensor and the approximate boost Killing vector field ^" 
as 



6Q 
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(1) 



Here 7 is the induced metric on the LCH and the dA 
integral is over a thin pencil of generators around the 
central generator connecting p to p^. 

We also need an entropy associated with the LCH. To 
begin with, taking a clue from the black-hole thermody- 
namics, it seems reasonable to guess that one should as- 
sociate the same expression of entropy to LCHs as that 
of the black holes in the theory. Let us then assume 
that the underlying theory of gravity obeys an equation 
Eab = 8 TT Tab and the entropy of any slice of LCH is same 
as that of the black-hole horizon in that theory given by. 
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(2) 



where p_b/4 is the entropy density of the black- hole hori- 
zon. For example, in general relativity ps would be equal 
to 1. Let us now consider how this entropy changes along 
the null congruence fc". The entropy change is 



SS = 
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Next, we expand the entropy change in a Taylor series 
around the terminal cross section at A = A/, 
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where the suffix A/ indicates that the quantity in the 
bracket is evaluated at the final cross section A = A; 



V- 



Since L^ dX X appearing in 5Q in Eq. (fl]) is of the same 

order in A/ asL ^ dX (A — A/), the second term of expan- 
sion in Eq. Q, for the validity of the Clausius relation 
we must have that the first term of expansion in Eq. (Hi) 
vanishes. By our construction of LCH, 9 is zero at A/ on 
the central generator. Hence a necessary condition for 
the validity of the Clausius relation for the LCH is that 



fdpB\ ^ 
\dXj,^. 



(5) 



In general relativity pB is 1 and Eq. ([5]) is trivially satis- 
fied. Hence, for general relativity we get from Eq. Q, 
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where in the second step we have used the Raychaudhari 
equation and the fact that 0{Xf) = = cr{Xf). Now 
using the Einstein equation and doing the A integral we 
see that Eq. ^ is equal to 5Q/T, where 6Q is as in 
Eq. (fTl) and T is the local Unruh Temperature, T = l/2n. 
Hence it is easy to prove the Clausius relation for LCHs 
in general relativity. For any other theory of gravity, 
however, the situation is different. For example, consider 
a theory of gravity described by the Lagrangian which is 
a function of the Ricci scalar f{R). The Wald entropy 
density associated with the black holes of such a theory 
is ^ f'{R), and the validity of Eq. ([S]) requires that 



dR 
dA 



= 0. 



But the spacetime is completely arbitrary and in general 
we can not set dR/dX = on the final cross section of 
the LCH. 

It is worth pausing here to recall the setting of the 
physical-process version of the first law of black-hole ther- 
modynamics: a black hole in an initial stationary state 
is perturbed by a small inflow of energy, finally it set- 
tles down to another stationary state and in the process 
there is some net change in its entropy |4l[T5l[T6j. Hence, 
for the physical processes involving black holes, the fi- 
nal cross-section is a part of a stationary spacetime, and 



therefore, all changes of the dynamical fields vanish in the 
asymptotic future. In particular, we can set dR/dX = 
on the final cross-section. 

This is the crucial geometric difference between the 
case of black holes and that of the LCHs. The future 
boundary condition for black holes is physically moti- 
vated from the stability and cosmic censorship hypothe- 
sis [17 . There is no compelling reason to impose such a 
condition on the LCHs. This is the difficulty in imposing 
Eq. (Isl) necessary for deriving the Clausius relation for 
the LCHs. 

It was suggested in Ref. [H] that one could impose a 
new equilibrium condition on the final cross-section of 
the LCH. Instead of requiring that 6 = Q a.t p one could 
choose a patch E such that the first term of the expansion 
in Eq. (|4| vanishes identically at p so that one doesn't 
need to impose Eq. ([5]) in addition to imposing 9 = 
at p. Then one could derive a non-equilibrium version of 
the Clausius relation that contains an extra term which is 
interpreted as the entropy-production term. In Ref. [S], 
this entropy production term is shown to be a separate 
contribution to the heat flux coming from the additional 
scalar degree of freedom present in f{R) gravity. Al- 
though such an interpretation offers an understanding of 
the entropy production terms for the specific case of f{R) 
gravity, it is not clear how to generalize it for a broader 
class of theories. In fact, it is explained in Ref. [13, why 
such an approach may not work beyond the simple case 
oi f{R) gravity 

A new proposal for the entropy density. — In this paper 
we follow a different approach. We believe that there is 
no compelling reason why the expression for the entropy 
density of stationary black holes should agree with that 
of the LCHs. We therefore seek a new candidate for the 
entropy associated with the LCHs. But this expression 
for the entropy density can not be completely arbitrary. 
There is some minimal set of properties one expects from 
this entropy density: it should of course satisfy the Clau- 
sius relation; as in the black-hole case it is desirable to 
have 9 = and aab = as the definition of the equi- 
librium slice of the LCH; for a stationary black hole the 
proposed entropy should agree with the entropy which 
satisfies the physical-process version of the first law for 
black holes. For example, for general relativity it should 
be proportional to the area of the horizon, for the Gauss- 
Bonnet gravity it should agree with the Jacobson-Myers 
expression |18) . 

With these guidelines in mind, our proposal is that the 
entropy density p of any slice E^ of LCH, located between 
p and po at an affine parameter A € [0, A/], is given by 



p = PB + {Xf - A) 



dps 
dX 
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psi^afefc^fc" - ^ - i^a^rfc^ ) , (7) 



where Pb/^ represents the expression for the entropy den- 
sity of the stationary black holes of the theory in ques- 



tion. The entropy of the shce is given by 



1 



dA^p. 



(8) 



Note that the entropy density of the final equilibrium 
slice is equal to pB , the entropy of the stationary black 
hole in the theory. 

Let us now check if our expression of the entropy den- 
sity leads to the correct thermodynamics of the LCH. The 
change in entropy as we evolve the system from the ini- 
tial slice at A = to the final equilibrium slice at A = Ay- 
is given by Eq. Q, with ps now replaced by p in Eq. ([7|, 
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dp 
dX 



(9) 



Now, noticing that ^(Ay) = because of our equilibrium 
condition, the first term in the big square brackets just 
gives {dp/dX)^ which is easily seen to be equal to from 
Eq. Q. The coefficient of (A— Ay) can be calculated in a 
straightforward fashion using the Raychaudhari equation 
and it gives simply {—Eabk°'k^)\f In the limit that the 
initial slice gets very close to the final equilibrium slice 
we get. 
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In the same limit we get, from Eq. (IT]) for the total heat 
flux through the null surface enclosed between A = and 

A = Ay, 



A/ 



dA^Tabk'^k' 



(11) 



Now using the local boost temperature T — 1/2it and 
the equation of motion Eab — SnTab, we see from equa- 
tions ( 10 1 and ( 11 ) that the Clausius relation SS = 5Q/T 
is satisfied at the LCH. 

Discussion. — Few comments on our proposal for the 
entropy density associated with the LCHs are in order 
now. First of all, in the case of general relativity, p^ = 1 
and Eabk°'k^ = Rabk°'k^ ■ Therefore, our entropy for- 
mula gives the usual area proportionality. Also, as we 
have already mentioned, on the final equilibrium slice 
our entropy density p agrees with pB , the entropy density 
of the stationary black-hole horizon in the theory. One 
might think that the previous statement is a tautology 
because on the final equilibrium slice one has A = Ay, 
so from Eq. ([7| the only surviving contribution to p is 



Pb- However, what is non-trivial is that for stationary 
black holes this should be true even for an arbitrary A, 
provided that the expression represented by Pb satisfies 
the physical-process version of the first law of black-hole 
thermodynamics. In order to see this, we note that all 
the derivatives of pB with respect to A vanish on any sta- 
tionary slice of the black-hole horizon. Furthermore, by 
stationarity 9, a and d9/dX on a stationary slice are all 
zero too, so Rabk°'k'' is zero by the Raychaudhari equa- 
tion. Finally, if pB satisfies a physical process law, there 
can not be any flux of matter stress tensor across the sta- 
tionary slice and hence Tabk"'k'' is zero which implies, by 
the equation of motion, that Eabk°'k'' is zero. Therefore, 
we see from Eq. ^ that for any stationary slice of the 
black- hole horizon, we have p = pb- This immediately 
shows that even for an arbitrary cross-section of a non- 
stationary black hole, our entropy formula gives p = pb 
as long as the perturbed black hole flnally settles down 
to a stationary configuration in the asymptotic future. 

We emphasize that it is not an aesthetic property that 
we demand from a proposed entropy density of the LCHs 
but it is essential for the consistency of the proposal with 
the black-hole thermodynamics. The reason is simple: 
the point p could be on any slice of the final stationary 
horizon of a black-hole and p would still be required to 
reduce to ps on this slice. From Eq. ^ this would be 
possible only if both the first-order and the second-order 
coefficient of (Ay — A) were to be zero for a stationary 
black-hole horizon, which is guaranteed to happen if pB 
satisfies the physical-process version of the first law for 
black holes. One could question the dependence of our 
p on the affine-parameter A and on our choice of the 
null vector /c". We should remember, however, in hy- 
drodynamics there is no derivation of the expression for 
a non-equilibrium entropy. One introduces the viscos- 
ity coefficients by hand to make sure that the entropy 
production is positive in the system's approach to the 
equilibrium. The same is true here. We have introduced 
terms by hand which ensure that the Clausius relation is 
satisfied. Note that these terms are independent of the 
additive and the multiplicative ambiguity in the choice of 
the affine parameter A. Furthermore, if the matter satis- 
fies the null energy condition then validity of the Clausius 
relation immediately implies the validity of the second 
law of thermodynamics for quasi-stationary approach to 
equilibrium. Therefore, in our case entropy production 
is guaranteed to be positive if the matter satisfies the 
null energy condition. It would be nice to understand 
if the A-dependent terms in our entropy density have an 
interpretation in terms of the hydrodynamic coefficients. 

We should compare our entropy density to that pro- 
posed in Ref . [13j . The situation in our proposal is better 
than the "Noetheresque" proposal of Ref. [13_. There the 
entropy density depended upon the choice of the approx- 
imate boost Killing vector field ^° and one had to make 
sure that the Killing identity is satisfied to an appropri- 
ate order for the derivation to go through. In contrast, 
our entropy density depends upon the null generator A;" 



of the LCH (relevant portion of) S which we chose to be- 
gin with such that its expansion and shear vanish at the 
equihbrium point p and we do not have to worry about 
the Kilhng identity at alL In both the cases though, the 
Clausius relation applies only to the complete patch of 
LCH between A = and A/. This is precisely analogous 
to what happens in the physical process version of the 
first-law for black holes. 

How sensitive is our derivation of the Clausius rela- 
tion to the choice of the equilibrium condition 6(Xf) = 0, 
a'^^{Xf) — 0? While d{Xf) = is essential for our deriva- 
tion, we can easily relax the condition <j(A/) — 0. This 
results in an extra term in Clausius relation that has a 
natural interpretation as an internal entropy production 
term ^ cr°'^aab whose coefficient can be interpreted as the 
shear viscosity [5]. For black holes in general relativity 
and f(R) gravity these are precisely the terms which ac- 
count for the tidal heating [9^ . To the best of our knowl- 
edge the expression for the shear viscosity in general the- 
ories of gravity is not known. 



Finally, we should mention that a physical interpre- 
tation of p has not emerged from our construction. We 
have a quantity that satisfies the Clausius relation for 
local causal horizons in all diffeomorphism invariant the- 
ories of gravity. But is there a derivation of this expres- 
sion from some basic principles? As far as this question 
is concerned, the situation is better in the proposal of 
Ref. [13] because that construction is based upon Wald's 
Noether charge entropy of the black holes in the theory. 
It remains to be seen if a similar interpretation is possible 
for our proposal of the entropy density of the LCHs. 
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